Non-Markovian dynamics of clusters during nucleation 
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Most theories of homogeneous nucleation are based on a Fokker-Planck-like description of the behavior of the 
mass of clusters. Here we will show that these approaches are incomplete for a large class of nucleating systems, 
as they assume the effective dynamics of the clusters to be Markovian, i.e., memoryless. We characterize these 
non-Markovian dynamics and show how this influences the dynamics of clusters during nucleation. Our results 
are validated by simulations of a three-dimensional Ising model with locally conserved magnetization. 
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Nucleation is the process where a stable nucleus spon- 
taneously emerges in a metastable environment. Excellent 
books and reviews exist on this topic; in his recent book, 
Kashchiev HI lists about 30 books and 40 review articles on 
nucleation. Examples of nucleation abound, for instance the 
formation of droplets in undercooled gasses and of crystals in 
undercooled liquids. The process is thermally activated and 
is key to understanding various subjects in biophysics, poly- 
mer physics, and chemistry. The physics behind it has long 
been studied and the simplest version is known as classical 
nucleation theory (CNT) 

In CNT the variations in the mass of a nucleus are described 
as a Markovian stochastic process in which single units attach 
and detach from the nucleus. The probability p(m, t) that a 
nucleus has mass m at time t evolves via the Fokker-Planck 
equation 00] 
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with R(m) the rate at which clusters of mass m grow to clus- 
ters of mass m + 1, F(m) the free energy of a cluster of mass 
77i and (3 the inverse temperature. Phenomenological expres- 
sions for R{m) and F(m) are provided to complete the the- 
ory. F(m) increases for small m, but decreases for large m, 
and it retains its maximal value at the so-called critical nu- 
cleus mass m c . Therefore, clusters with mass below m c tend 
to shrink, whereas clusters larger than m c tend to grow. Start- 
ing at time t = with cluster mass m(0) = m c , the mass 
evolves diffusively in time: its mean square change, averaged 
over all realizations, (Am(t) 2 ) = ((m(t) — m c ) 2 ), grows lin- 
early with time for small deviations from the critical mass. 

Although, qualitatively, there is ample experimental verifi- 
cation of the predictions of CNT, quantitatively, the nucleation 
rates it predicts may differ from experiment by 5 to 10 orders 
of magnitude [5] . CNT works really well if the nucleation pro- 
cess is correctly described by a Markov process with units at- 
taching and detaching as uncorrelated events. This is more or 
less the case in, for example, the Ising model with spin flip dy- 
namics. However, in the presence of a local conservation law, 
which, for example, is the case in binary mixtures of fluids or 
gasses and the Ising model with local spin exchange dynam- 





FIG. 1 : Illustration of memory effects in the Ising model with local 
spin exchange dynamics (exaggerated). Although they have equal 
masses, the cluster in the left panel is likely to grow as it has just 
shrunk, hence has higher density surrounding it. In the right panel 
this is vice versa. 



ics, CNT shows serious shortcomings due to the neglect of 
memory effects by describing the cluster growth as a Markov 
process. The most dominant of these usually is the strong cor- 
relation in time between attachment and detachment events: 
after a particle detaches from the nucleus, it remains in its 
neighborhood for a relatively long time and is therefore likely 
to reattach. In fact, the theory of Brownian motion [6] states 
that in two dimensions, every detached particle will eventually 
return to the cluster, possibly after an extremely long time, 
and that in three dimensions, a fraction of the detached parti- 
cles will never return. Figure [T]illustrates this memory effect 
in the Ising model with local spin exchange dynamics. Two 
clusters of equal mass are shown. The cluster in the left panel 
has just shrunk, consequently has a higher density of particles 
surrounding it, which enhances the probability that the cluster 
will grow; in the right panel, the opposite is happening and the 
cluster is more likely to shrink. These memory effects result in 
behavior quite different from CNT. The relation between the 
driving force —OF /dm and the drift velocity is more com- 
plicated and the diffusion is anomalous, i.e. the mean square 
displacement scales non-linearly with time. 

We develop a theory for the variation in the cluster mass, 
where these memory effects are taken into account and show 
that on different time scales different types of dynamical 
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(<J£r(t)*£r(f )) = PrBjJf (m(t))*(t - f)> 



(4) 



FIG. 2: The four effects causing variations in the cluster mass: 1) 
the emission of returning particles, 2) the emission the non-returning 
particles, 3) the absorption of non-returning particles, and 4) the ab- 
sorption of returning particles. 



behavior take place. We concentrate on states in three- 
dimensional systems with a small gradient in free energy, 
which is the case for, for instance, near-critical nuclei. In that 
case, three distinct regimes show up. On very short time scales 
the mean square growth of a cluster scales linearly with time 
like in an ordinary diffusive process. This regime is followed 
by a regime of anomalous diffusion with exponent of one half 
(i.e. (Am(t) 2 ) ~ t 1 / 2 ) and finally on large time scales linear 
growth occurs again, but with a much smaller prefactor than in 
the first regime. We validate our theory with extensive simula- 
tion results on the mean square growth of near-critical nuclei 
in a three-dimensional Ising model with local spin exchange 
dynamics. 

In our model, we consider a cluster in a dilute environment 
(like in figureQ]). The variations in the cluster mass are caused 
by four effects, namely 1) the emission of particles from the 
cluster, that are returning (with return probability p r ), 2) the 
emission of non-returning particles, 3) the absorption of par- 
ticles from far away (i.e. not returning from previous emis- 
sion), and 4) the absorption of particles that are returning. 
These four effects are illustrated in figure[2] The first three ef- 
fects are independent and are described by random functions 
£ r (t), £,n r (t) and £n r (t), respectively. Each of these random 
functions consists of a series of delta functions at Poisson dis- 
tributed random times, so that the cluster mass m(t) is integer 
at all times. Their average values obey 



(&•(*)> -<6r(*)> =«("*(*)), 



(2) 



with v(m) the systematic growth of a cluster of mass m. The 
variations 6£ a (t) = £ a (t) — (£ Q (i)) obey 

X>&(*)*e£(*0) = 2 (1 -Pr)D M (m(t))S(t - f), (3) 



with Djvf(m) the diffusion coefficient describing the short 
time mass variations of a cluster of mass m. Equations (O- 
dUi fully determine the three random functions 

The absorption of returning particles is correlated to their 
emission: a returning particle emitted at time r is returning at 
time r + T T , with T T described by a return time distribution 
/j. Putting this together, one may describe the evolution of a 
cluster by the stochastic differential equation: 



™(*)=£U*)-Cr(*)-£r(i) 
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drd(t 



(5) 



The first two terms are basically CNT and the addition of the 
other two terms is new. This stochastic differential equation is 
too hard to be solved in general and therefore we restrict our- 
selves first to near-critical clusters. On average critical clus- 
ters absorb equally many particles as they emit, so we assume 
v(m) = 0. Furthermore, we assume that the diffusion coeffi- 
cient is m-independent: Dm(iti) = Dm- The average growth 
(Am(t)) is then vanishing and the mean square growth can be 
calculated from 



(Am(() 2 ) = dr dr' (m(r)TO(r')) , 



(6) 



which, after substituting equation §5^ twice and using (O and 
©, results in the following expression 



(Am(() 2 } 



2D M [t-Pr \dr{t - t)h{t) 



(7) 



The first term in this equation is ordinary diffusion, which is 
suppressed by the last term, since emission and absorption are 
correlated at larger time scales. 

Next we investigate various limits of the equation. We 
assume that the time scale of spontaneous fluctuations (i.e. 
t ps 1/ Dm) is much smaller than the time scale on which par- 
ticles are typically returning. For small times the leading term 
in equation results in 



(Am(i) 2 ) fa 2D M t 



(8) 



At these time scales attachments and detachments occur at dif- 
ferent places on the cluster's surface as independent events, 
and therefore the mean square cluster growth scales linearly 
with time. 

To investigate (Am(() 2 ) at larger times we have to specify 
fi(t) in more detail. We assume that a detached particle makes 
a three-dimensional random walk with diffusion coefficient 
Db- Furthermore, it starts at a distance 5R separated from 
the cluster, which we consider to be a sphere of radius R. This 
gives for the return probability p r = R + SR and for the return 
time distribution 101 



SR exp (-5R 2 /AD B t) 



(9) 
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Asymptotically at large times p(i) ~ t~ 3 / 2 , but it is cut off at 
small times by the exponential. Using this distribution to cal- 
culate the mean square cluster growth for larger times results 
in 

(Am(t) 2 ) w 2(1 -p r )D M t + const y/t. (10) 

At large time scales, the first term dominates, and the clus- 
ter dynamics is determined by particles being emitted to and 
absorbed from far away with rate 2(1 — p r )Dj^. At those 
large times, the change in cluster size is determined by non- 
returning absorbed and emitted particles, and those particles 
can be treated as independent events. At intermediate time 
scales, however, the anomalous diffusion dominates if p r is 
close to one, which it typically is. This behavior is caused 
by emitted particles returning, much comparable to single file 
diffusion QSJ] . Note that if the return probability is one, the 
anomalous diffusion lasts forever. This may happen, for ex- 
ample, in many two-dimensional systems, although the long 
time behavior is different from equation ( [Tol l, due to a differ- 
ent return time distribution p(t). 

To validate the presented theoretical picture, simulations 
were performed on a three-dimensional Ising model with lo- 
cal spin exchange (Kawasaki) dynamics with Metropolis ac- 
ceptance probabilities J9l LlOfl , a prototypic system to study 
nucleation. The Hamiltonian is given by 

H = -J^SiSj + h^Si, (11) 

(id) » 

with the first sum over all pairs of adjacent sites. A three- 
dimensional cubic lattice of size 32 x 32 x 32 with periodic 
boundary conditions has been used for the simulations. The 
temperature is chosen as ksT = 2.5 J, well below the critical 
temperature of fcgT c ss 4. 5 J [10]. An oversaturated initial 
configuration is first brought into equilibrium under constant 
magnetization, resulting in a single large cluster of m c spins 
in equilibrium with its surroundings with a density of 1.2 per- 
cent. Next, besides the spin exchange moves, we also perform 
spin flips in three strips of the box (one in each principal direc- 
tion), as far away from the cluster as possible. These spin flip 
moves mimic an infinite reservoir of up-pointing spins, and 
render the cluster instable. The strength of the external field 
h is fixed at the value giving the initial cluster a critical size, 
hence zero average growth rate. The free energy as a function 
of cluster mass is plotted in figure [3] before (left) and after 
(right) lifting the constraint on the magnetization. 

The time evolution of this cluster is then measured and the 
resulting mean square change in cluster mass as a function of 
time, averaged over about 70,000 realizations, is shown in fig- 
ure |4j For comparison, our theoretical estimate of (Am(i) 2 ) 
is plotted with the simulation data. Our theory requires as in- 
put the function p[t), with its parameters Db, R and SR, and 
the parameters Dm and p r . They are chosen as follows. The 
distribution p(t) is as in equation Q, multiplied by an expo- 
nential, since the system has a finite size and long return times 
are exponentially surpressed. The result depends very little on 
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FIG. 3: Free energy landscapes of a finite-sized Ising model, as a 
function of cluster mass; in the left panel with conserved magnetiza- 
tion and in the right panel without. The density in the left panel and 
the external magnetization in the right panel are tuned, such that the 
extrema of the potentials coincide. 



the chosen exponent. Furthermore, Db = |, in accordance 
with the definition of time in our model, and R is chosen such 
that m c = ^irR 3 . The parameter 5R is the initial distance be- 
tween the cluster and a spin which has just detached from it. If 
the spin detaches radially, this distance equals the lattice spac- 
ing, but the effective distance can be smaller by a factor of up 
to two in other directions. We chose SR = , since that is 
the average distance between a spherical cluster and the sites 
neighboring it on the square lattice. The diffusion coefficient 
Dm in equation (0 and subsequently the return probability p r 
are fitted. 

The theory captures the general trend well. For short times 
ordinary diffusion is observed in the simulation data, after 
which anomalous diffusion with exponent one half is clearly 
present, as indicated by the dashed line in figure [4] For large 
times, however, the cluster is so far out of equilibrium that the 
gradient in the free energy is non-zero, so that (Am(t)) ^ 
and (Am(t) 2 ) grows super-linearly with time. We verified 
that the same results hold, qualitatively, for different values of 
the temperature. 

Next, we turn to the consequences of these memory effects 
on nucleation times. To obtain nucleation times, equation 
(0 should be solved in the presence of a free energy barrier, 
which results in a systematic growth v(m) in equation ©. 
The phenomenological expression for the free energy from 
CNT J2, or quadratic approximation could be used. This 
work is still in progress and therefore we resort to scaling ar- 
guments in the meantime. 

Instead of the nucleation time, we focus on the vaporiza- 
tion time (i.e. the time it takes for a critical nucleus to vapor- 
ize); these two times are connected via detailed balance by 
Tnuci ~ e l3AF T vav , hence T vap is the pre-exponential factor of 
the nucleation time. This vaporization time is mainly dictated 
by the time a critical cluster resides in a region near the top 
of the free energy barrier; after that it vaporizes relatively fast 
due to the gradient in the free energy. We define this region 
near the top as the region where the free energy is less than 
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krjT below the maximum value F{m c ) and call its width 8m. 
In CNT the residence time for this region, often referred to as 



the Zeldovich factor [11], scales as T res ~ 8m 2 . If, in equa- 



tion ( [Tol l, the anomalous diffusion is taken over by the normal 
diffusion at the residence time, the memory effects only re- 
sult in rescaling the diffusion coefficient by a factor of 1 — p r . 
However, if the cluster dynamics show sub-diffusive behavior 
up to this residence time, then (Am(i) 2 ) ~ \fi, so that the 
residence time is proportional to 8m 4 . Which behavior occurs 
depends on numerous variables, such as the shape of the free 
energy landscape and the mobility of detached particles. 

Another view on the consequences of memory effects on 
the nucleation time is that, compared to CNT with rescaled 
diffusion coefficient Z? e ff = (1— p r )Dj^, anomalous diffusion 
is present in addition to CNT's ordinary diffusion. This results 
in increased mass fluctuations around the critical mass, hence 
decreases the residence time and therefore also the nucleation 
time. 

In summary, we have demonstrated that memory effects 
are playing an important role in the dynamics of nucleation. 
The time evolution of nucleating clusters shows consequently 
anomalous diffusion. These memory effects can be captured 
by a simple stochastic differential equation (0, which gives 
measurable results for critical clusters. These effects proba- 
bly also have impact on the nucleation times. Analytic treat- 
ment of this equation in the presence of a free energy barrier 
is therefore an interesting topic for future research and might 
lead to new quantitative predictions for nucleation times. 

We thank Henk van Beijeren for useful discussion. 
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FIG. 4: Measurements of (Am(t) 2 )/t in the Ising model, plotted 
together with the theoretical value of equation 10 and a linear fit with 
a slope of minus one half. A zero slope indicates ordinary diffusion. 



